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Abstract 

We apply a data-driven method to STAR Au+Au collisions at V*aw = 200 GeV to iso- 
late A^-dependent and A^-independent correlations by using two- and four-particle Q-cumulant 
v„ measurements. The A^-independent part, dominated by flow, is found to be 77-independent 
within the STAR TPC of +1 unit of pseudo-rapidity. The A^-dependent part may be associated 
to nonflow, and is responsible for the A77 drop in the measured two-particle V2 cumulant. We 
combine our result to four- and six-particle cumulants to gain further insights on the nature of 
flow fluctuations. 



1. Introduction 

Heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) provide a means to study 
the Quark Gluon Plasma (QGP). In a non-central collision, the transverse overlap region is 
anisotropic. The energy density gradient converts the initial coordinate-space anisotropy into 
the final momentum-space one, generally called flow. As the system expands, the coordinate- 
space anisotropy diminishes. Hence, flow is sensitive mostly to the early stage of the collision 
[1 1. Through measurements of anisotropic flow and comparison to hydrodynamic calculations, 
properties of the early stage of the collision system may be extracted. For example, the shear 
viscosity to entropy density ratio (77/s) of the QGP was found to be not much larger than the 
conjectured quantum limit of 1 /4n 0. 

The momentum-space anisotropic flow can be characterized by the Fourier coefficients of 
the particle azimuthal distribution relative to a symmetry plane (participant plane) Q. However, 
the participant plane is not experimentally accessible. As a proxy, the event plane is constructed 
from particle momentum. Similarly, two- and multi -particle correlations are also used to quantify 
anisotropy. These various flow measurements are therefore contaminated by intrinsic particle 
correlations unrelated to the participant plane, generally called nonflow J4|. Although nonflow 
at low transverse momentum (p T ) may be small, it was found that 77/ s is very sensitive to flow 
0. Therefore, it is important to separate the nonflow contributions from flow measurements. 
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2. Measurement Methods 



This analysis attempts to separate flow and nonflow in a data-driven manner [6|. The method 
utilizes two-particle cumulant with one particle at r\ a and another at r\p , and four-particle cumu- 
lant with two at r\ a and two at r\p. We have 

V{2) = «2» = v(r, a )v(r,p) + cr(r, a )cr( m ) + cr'iAjj) + 8(Arj), (1) 

V{4] = V2«2» 2 - «4» * vCi/Xij,) - crinaMrjp) - (r'(Arf), (2) 

where Ar\ - r\p- rj a , and we require j] a < Tjp < 0. The average flow (v) is only a function of 77. 
cr and <t' are the flow fluctuation, which may have A/7-dependence component. Nonflow (6) is 
generally a function of A77. 

We take the difference between cumulants V„{2} at (T] a , rjp) and {rj a , —rjg). We do the same 
for V„{4). For symmetric collision systems, the A^-independent terms in Eq. ([T} and (|2]i cancel 
out, leading to Eq. |3]l and Q: 

AV{2) = V{2}(rj a , rjp) - V{2}(77 ff , -r, p ) = y{2}(A/ 7l ) - V{2}(A^ 2 ) = At/ + M, (3) 
AV{4) = V{4}(;7 ff , ^) - V{4}(7, ff , -77^) = y{4}(A/ 7l ) - V{4}(At7 2 ) * -Act', (4) 

where Act' = <x'(At/i) - ct'(At]2),A6 = 5(A^i) - 5(A^), and A/71 = rjp - rj a , A772 = -77^ - 770-. 
Assuming the same flow fluctuations cr' in both V{2} and V{4), we can extract information of 
nonflow (6) and A^-dependent flow fluctuation (cr'). Using Eq. ([Tji and Q, we can further obtain 
v and cr. 

3. Data Analysis and Results 

We analyze the data in Au+Au collisions at V*mv = 200 GeV from year 2004 measured 
by the STAR Time Projection Chamber (TPC) within \rj\ < 1.0. The Q-cumulant method QE) 
is used with the acceptance correction in calculation. Azimuthal moments are weighted by unit 
(not pair/quadrant multiplicities). 

Figure [TJ (a), (b) and (c) show the cumulant results in Eq. ([T]) and Q as a function of 77,,., 77^, 
and (d), (e) and (f) show the quantities in Eq. ([3]) and (|4]) for 20-30% centrality. AV2{4) is shown 
in Fig.[TJ(d), and is consistent with zero, suggesting that elliptic flow fluctuations are independent 
of A77 within the TPC acceptance. This is used in subsequent analysis of data. We further assume 
the triangular flow fluctuations are also independent of A77 within the TPC acceptance. With this 
taken into account, AV{2) will only have nonflow component. AV2{2) and AV3{2) is shown in 
Fig.[T](e) and (f), respectively. 

As the dashed lines shown in Fig. [T] (e) and (f), the nonflow difference between A771 and 
A7/2 appears linear in A772 - A/71 , with a slope that is independent of A7/1. The intercept of the 
linear dependence decreases quickly with A/71, suggesting a short-range exponential component 
in nonflow. In other words, the nonflow difference can be described by 

6(Atji) - 6(Atj 2 ) = a[e~ - e~^] + A[e ^ 2 -e^ 2 ]. (5) 

A single fit of Eq. <|5j with four free parameters a, b, A, cr is applied to all the data points, yielding 
the solid curves in Fig.[TJ(e) and (f). The fit results give the following nonflow parameterization: 

-A;; -A?T 

6(Aj]) - ae i' + Ae 2^ . (6) 
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Figure 1 : (a) The second harmonic four-particle cumulant. (b,c) The second and third harmonic 
two-particle cumulants. (d) Difference in V2{4) between (j] a , rjp) and (i] a , -rip) as a function of 
A772 - A/71 = \2rjij\. (e,f) Difference in V2{2) and V^{2} for given A/71 values, respectively. Data 
are from the 20-30% centrality of Au+Au collision at -\Js^^ — 200 GeV. 

Subtracting Eq. |6) from the raw two-particle cumulant gives the decomposed 'flow' results in 
Fig. [2] (a) and (b). Here the decomposed flow includes the flow fluctuation effect (v(rj a )v(rij3)) = 

Figure|2](a) and (b) shows that vv + crcr is independent of 77 for both \'i and V3, and Fig.[T](a) 
shows V2V2-0"2CT2 is independent of r\. This suggests that v>2 and its fluctuation are independent of 
77. The projections of Fig.[2](a), Fig.[T](a) and (b) are shown in Fig.|2](c). The difference between 
the raw two-particle cumulant V2{2) and the decomposed 'flow' (v 2 ,) is the nonflow effect. We 
find that 62lv\ is about 4%, and a\/v\ is about 13% for the centrality 20-30%. Figure [3] shows 
the decomposed 'flow' and nonflow as a function of centrality in Au+Au V*aw - 200 GeV. 

The flow fluctuation may also be assessed by four- and six-particle cumulants. Assuming a 
Gaussian fluctuation of the flow magnitude, we have v{6}/v{4) * 1 - cr 6 /3v 6 , when cr 2 /v 2 « 1 
[8|. On the other hand, independent Gaussian fluctuations of v x ,v y would give v{6) — v{4) [9|. 
We show in Fig. |4] V2{6}/v2{4) vs the fluctuation results from our decomposition method. The 
systematic errors are estimated by the different moment weight methods and non-uniform ac- 
ceptance correction. The v'2 results, while favoring Gaussian fluctuations in the flow magnitude, 
may be also consistent with Gaussian fluctuations in the individual flow components. A firmer 
conclusion on the nature of flow fluctuations await for more statistics. 

4. Summary 

We analyze two- and four-particle Q-cumulants between pseudo-rapidity bins in Au+Au 
collisions at 200 GeV from STAR. Exploiting the collision symmetry about mid-rapidity, we 
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Figure 2: (a) {v\}. (b) {v\) as a function of (T] a , 77^). The parameterized nonfiow has been sub- 
tracted, (c) the projection of 77 dependence of the decomposed flow {v\) and the raw two- and 
four-particle cumulants of V2 for the 20-30% centrality of Au+Au collision at V*aw = 200 GeV. 
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Figure 3: the decomposed flow and nonflow Figure 4: V2{6}/v2{4} and cr\l\\. The curves 
with 77-gap > 0.7 as a function of centralities are model predictions of ref. |9 1 and J8) 



isolate the A/7-dependent and A^-independent correlations in the data. Elliptic flow fluctuations 
appear independent of A77 within the TPC acceptance, as deduced from the four-particle cumulant 
measurement. The method does not make assumptions about the rj dependence of flow. Our 
isolated A^-independent part of our data is dominated by flow, with a relatively small away-side 
nonflow, and is found to be 77-independent within the STAR TPC of ± 1 unit of pseudo-rapidity. 
The A77-dependent part is nonflow and is found to decrease with A77. In the 20-30% Au+Au 
collisions at \rj\ < 1, and pj < 2 GeV/c, the elliptic nonflow relative to flow ^/(vf) is found to be 
4% and the relative fluctuation <r|/(v^) is 13%. The flow fluctuation in conjunction with the 6- 
and 4-particle cumulants flow ratio seems to favor Gaussian fluctuations in the flow magnitude. 
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